Non-standard Neutrino Oscillations at Icecube 
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Abstract 

In this talk I review the potential of Icecube for revealing physics beyond the standard 
model in the oscillation of atmospheric neutrinos [1]. 



1. Introduction 



J> | With its high statistics data [2] Super-Kamiokande (SK) established beyond 

doubt that the observed deficit in the ^t-like atmospheric events is due to oscillations, 
a result supported by the K2K and MINOS long-baseline (LBL) experiments [3,4]. 

Mass oscillations are not the only possible mechanism for atmospheric (ATM) 
Vp, — * v T flavour transitions. These can be also generated by a variety of nonstan- 
\& ■ dard physics characterized by the presence of an unconventional v interaction that 

mixes neutrino flavours [5] . Examples include violations of the equivalence principle 
p ■ (VEP) [6,7], non-standard neutrino interactions with matter [8], neutrino couplings 

p { to space-time torsion fields [9], violations of Lorentz invariance (VLI) [10] and of 

CPT symmetry [11,12]. In contrast to the E energy dependence of the conventional 
| oscillation length, new physics can produce neutrino oscillations with wavelengths 

p . that are constant or decrease with energy. [13,14]. At present these scenarios can- 

not explain the data[15] and a combined analysis of the ATM and LBL data can 
be performed to constraint them even as subdominant oscillation effects [16]. 
IceCube, with energy reach in the 0.1 ~ 10 4 TeV range for ATM neutrinos, is 



the ideal experiment to search for new physics. For most of this energy interval 
standard Am 2 oscillations are suppressed and therefore the observation of an an- 
gular distortion of the ATM neutrino flux or its energy dependence provide a clear 
signature for the presence of new physics mixing neutrino flavours [1] . 

2. Propagation in Matter of High Energy Oscillating Neutrinos 

We concentrate on v il -v T flavour mixing mechanisms for which the propagation 
of u's (+) and P's (— ) is governed by the following Hamiltonian [12]: 



Am 2 is the mass-squared difference between the two neutrino mass eigenstates, er„ 
accounts for a possible relative sign of the new physics (NP) effects between i/s and 
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z/s and A<5„ parametrizes the size of the NP terms. By r\ n we denote the possible 
non-vanishing relative phases. 

cos6> sin0\ / cos£„ sin£„e ±M '™ \ 

U 9 = | , V Ut±Tln = ?n ; (2) 

— sin 9 cos 9 1 \ — sm£ n e TlVn cos£„ / 

If NP strength is constant along the neutrino trajectory the oscillation probabil- 
ities take the form [12]: 
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K = y/l + R\ + 2R n (cos 20 cos 2£„ + sin 26* sin 2£ n cr] n ) , i?„ = a. 
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where, for simplicity, we have assumed scenarios with one NP source characterized 
by a unique AS n . crj n — cos r\ n 

Eq. (1) describes, for example, flavour mixing due to new tensor-like interactions 
for which n = 1 leading to a contribution to the oscillation wavelength inversely 
proportional to the neutrino energy. This is the case for i^'s and v T 's of different 
masses in the presence of violation of the equivalence principle due to non- universal 
coupling of the neutrinos, 71 7^ 72 to the gravitational potential cj> [6], so ASi = 
2\(f>\(ji — 72) = 2|0|A7. v\ and v 2 are related to and v T by a rotation £1 = t; vep . 

For constant potential <j> 1 this mechanism is phenomenologically equivalent to the 
breakdown of Lorentz invariance resulting from different asymptotic values of the 
velocity of the neutrinos, ci ^ c 2 , A<5i = (ci — c 2 ) = 5c/ c, with v\ and v 2 being 
related to and v T by a rotation £1 = £ v u [10]. 

For vector-like interactions, n = 0, the oscillation wavelength is energy-independent. 
This may arise, for instance, from a non-universal coupling of the neutrinos, k\ ^ k 2 
so A<5 = Q{ki—k 2 ) {v\ and v 2 is related to the and ^ T by a rotation £ = £q), to 
a space-time torsion field Q [9] . Violation of CPT resulting from Lorentz- violating 
effects such as the operator, ^l^I^L' a ^ so leads to an energy independent contri- 
bution to the oscillation wavelength [11,12] which is a function of the eigenvalues of 
the Lorentz violating CPT-odd operator, 6 i; A<5 = b\ — b 2 , and the rotation angle, 
Co = £(2pt, between the corresponding eigenstates Vi and the flavour states v a . 

For most of the neutrino energies considered here, Am 2 oscillations are suppressed 
and the NP effect is directly observed. Thus the results will be independent of the 
phase r] n and we can chose the NP parameters in the range A5 n > 0, < £„ < 7r/4. 

The Hamiltonian of Eq. (1) describes the coherent evolution of the v^-Vt ensem- 
ble for any neutrino energy. High-energy neutrinos propagating in the Earth can 
also interact inelastically with the Earth matter either by charged current (CC) and 
neutral current (NC) and as a consequence the neutrino flux is attenuated. This 
attenuation is qualitatively and quantitatively different for iVs an d z^'s. v^s are 
absorbed by CC interactions while v T 's are regenerated because they produce a r 
that decays into another tau neutrino before losing energy [17]. As a consequence, 
for each v T lost in CC interactions, another v T appears (degraded in energy) from 
the r decay and the Earth never becomes opaque to v' T s. Furthermore, a secondary 
flux of D^s is also generated in the leptonic decay r — > iiv^v T [18]. 
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Attenuation and regeneration effects of incoherent neutrino fluxes can be consis- 
tently described by a set of coupled partial integro-differential cascade equations 
(see for example [19] and references therein) or by a Monte Carlo simulation of the 
neutrino propagation in matter [17,18,26]. For astrophysical v's, because of the long 
distance traveled from the source, the oscillations average out and at the Earth the 
neutrinos can be treated as an incoherent superposition of mass eigenstates. 

For ATM neutrinos this is not the case because oscillation, attenuation, and 
regeneration effects occur simultaneously when the neutrino beam travels across 
the Earth's matter. For conventional neutrino oscillations this fact can be ignored 
because the neutrino energies covered by current experiments are low enough for 
attenuation and regeneration effects to be negligible. But for non-standard scenario 
oscillations, future experiments probe high-energy neutrinos for which the attenu- 
ation and regeneration effects have to be accounted for simultaneously. 

In order to do so it is convenient to use the density matrix formalism to describe 
neutrino flavour oscillations. The evolution of the neutrino ensemble is determined 
by the Liouville equation for the density matrix p(t) = v{t) <g> v{t)^ 



where H is given by Eq. (1). The survival probability is given by P^{t) = Tr\n. v p(t)], 
where IT^ = v^®v^ is the state projector, and with initial condition p(0) — IT^ . 
An equivalent equation can be written for the v density matrix. 

In this formalism attenuation effects due to CC and NC interactions can be 
introduced by relaxing the condition Tr(p) = 1. In this case 



where [Xf nt (E, t)}^ = [Ag c (£7, t^+^ciE, t)}~\ [Ag c (£7, t)}^ = n T (x) a% c (E), 



and [Anc(-E') t)]' 1 = n T (x) 0nc (E) (jit(i) is the number density of nucleons at the 
point x = ct). 

v T regeneration and neutrino energy degradation can be accounted for by coupling 
these equations to the shower equations for the r flux, F T (E T , t) (we denote by F 
the differential fluxes d(p/(dE d cos 9)). For convenience we define the neutrino flux 
density matrix F v {E 1 x) = F„ (E, Xo)p(E, x = ct) where F v (E, xq) is the initial 
neutrino flux: 



dp 



i[R,p], 



(3) 



dt 






dF T (E T ,t) 



1 



F T (E T ,x) 



dx 



Kee(Er,x) 
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X T dec {E T ,x) = 7 r cr T . r r is the r lifetime and j r = E T /m T is its gamma factor. 

dN NC (El,E„) = 1 da NC (El,E„) , dN cc {E v ,E T ) = 1 d^ c (E v ,E T ) , 
dE„ _ o-nc(B^) <i£ T — <rJ, c (_E„) dE T Uti 

easily computed. The r decay distributions dN ^^j,E u ~) an d dN doc {E T ,E v ) can ^ g 
found in Rcfs. [19,27]. 

The third term in Eq. (5) represents the neutrino regeneration by NC interactions 
and the fourth term represents the contribution from v T regeneration, v T — > t~ — > 
&V> describing the energy degradation in the process. The secondary f M flux from v T 
regeneration, v T — > r+ — > j/ r /i+ z^, is described by the last term where we denote 
by over-bar the energies and fluxes of the t + . Br^ = 0.18 is the branching ratio for 
this decay. In Eq. (6) the first term gives the loss of taus due to decay and the last 
term gives the r generation due to CC v T interactions. In writing these equations we 
have neglected the tau energy loss, which is only relevant at much higher energies. 

An equivalent set of equations can be written for the v flux density matrix and 
for the r + flux. Both sets are coupled due to the secondary v flux term. 

We solve this set of ten coupled evolution equations that describe propagation 
through the Earth numerically using the matter density profile of the Preliminary 
Reference Earth Model and obtain the neutrino fluxes in the vicinity of the detector 
%£§ = Tv[F v (E,L = 2Rco S e)IL a ] . 

In Fig. 1 we illustrate the interplay between the different terms in Eqs. (5) and (6). 
The figure covers the example of VLI-induced oscillations with 5c/ c = 10~ 27 and 
maximal £ v u mixing. The upper panels show the final and v T fluxes for vertically 
upgoing neutrinos after traveling the full length of the Earth for the initial condi- 
tions d^{vn)o/dE v = d${vn)a/dE v oc E^ 1 and d§(v T ) /dE v = d<&{v T )o/dE v = 0. 

The figure illustrates that the attenuation in the Earth suppresses the neutrino 
fluxes at higher energies. The effect of the attenuation in the absence of oscillations 
is given by the dotted thin line in the left panel. Even in the presence of oscillations 
this effect can be well described by an overall exponential suppression [27,24] both 
for v^'s and the oscillated v T 's. In other words, we closely reproduce the curve for 
"oscillation + attenuation" simply by multiplying the initial flux by the oscillation 
probability and an exponential damping factor: 

^°ifj'fr g) = dJ^h p ^ E > L = 2Rcos0 ^ eM-x(e)(aME)+^ cc (E))] , 

dEdcosv dEacostt 

(7) 

where X(9) is the column density of the Earth. 

The main effect of energy degradation by NC interactions (the third term in 
Eq. (5)) that is not accounted for in the approximation of Eq.(7) is the increase 
of the flux in the oscillation minima (the flux does not vanish in the minimum) 
because higher energy neutrinos end up with lower energy as a consequence of the 
NC interactions. The difference between the dash-dotted line and the dashed line is 
due to the interplay between the v T regeneration effect (fourth term in Eq. (5)) and 
the flavour oscillations. As a consequence of the first effect, we see in the right upper 
panel that the v T flux is enhanced because of the regeneration of higher energy ^ T 's, 
v T (E) — > t~ — > v T (E' < E), that originated from the oscillation of higher energies 
Vfjjs. In turn this excess of ^ T 's produces an excess of v^s after oscillation which 
is seen as the difference between the dashed curve and the dash-dotted curve in 
the left upper panel. Finally the secondary effect of v T regeneration (last term in 
Eq. (5)), v T {E) — > t+ — > [i + v T v^iE' < E), results into the larger flux (seen in 
the left upper panel as the difference between the dashed and the thick full lines). 
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This, in turn, gives an enhancement in the v T flux after oscillations as seen in the 
right upper panel. 




E v [GeV] E v [GeV] 



Fig. 1. Vertically upgoing neutrinos after traveling the full length of the Earth taking into account 
the effects due to VLI oscillations, attenuation in the Earth, v T regeneration and secondary v T 
regeneration (see text for details). 

The lower panels show the final fluxes for an atmospheric-like energy spectrum 
d^{v li )o/dE v = d^(D^) /dE v oc E~ 3 and d${v T ) /dE u = d<S>(D T ) /dE u = 0. In 
this case regeneration effects result in the degradation of the neutrino energy and 
the more steeply falling the neutrino energy spectrum, the smaller the contribution 
to the total flux. As a result the final fluxes can be relatively well described by the 
approximation in Eq.(7). 



3. Example of Physics Reach: VLI-induced Oscillations 

The expected number of v ^ induced events at IceCube can be obtained by a 
semianalytical calculation as: 

/l />OG />00 />00 />00 />OC 

dcosfl / dl' min / dl / dEl» / dEl / dE v (8) 
1 JO ~' ; min ^"V "' £ V n J E % 

dE^dcoa 6 ^ s ^ ne differential muon neutrino neutrino flux in the vicinity of the de- 
tector after evolution in the Earth matter obtained as described in the previous 
section. We use the neutrino fluxes from Honda [20] extrapolated to match at 
higher energies the fluxes from Volkova [21]. At high energy prompt neutrinos from 
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charm decay are important and it is evaluated for two different models of charm 
production: the recombination quark parton model (RQPM) developed by Bugaev 
et al [22] and the model of Thunman et al (TIG) [23] that predicts a smaller rate. 
d dE" (Ev>Eu) is the differential CC interaction cross section producing a muon of 

energy E®. T is the exposure time of the detector. Equivalently, muon events arise 
from interactions that are evaluated by an equation similar to Eq.(9). 

After production with energy E^, the muon ranges out in the rock and in the 
ice surrounding the detector and looses energy. We denote by F(E°,E^ n ,l) the 
function that describes the energy spectrum of the muons arriving at the detector. 
We compute the function F(E®, E^ n , I) by propagating the muons to the detector 
taking into account energy losses due to ionization, bremsstrahlung, e+e~ pair 
production and nuclear interactions according to Ref. [24]. 

The details of the detector are encoded in the effective area ^4°// f° r w hich we 
make a phenomenological parametrization to simulate the response of the IceCube 
detector after events that are not neutrinos have been rejected ( referred to as "level 
2" cuts in Ref. [25]). The explicit form of A° eff cn be found in Rcf.[l]. 

Together with ^-induced muon events, oscillations also generate /x events from 
the CC interactions of the v T flux which reaches the detector producing a r that 
subsequently decays as r — > ^VfiV T and produces a /x in the detector: 



/l r-OO r-OO r-OO r-OQ r-OO r-OC 

dcosO / dl' min / dl dE* n / dE° / dE T / dE v (9) 

-1 JO ' l 'min •' m ^ •' E ^ D 



where — d ° c dE o' " can be found in Ref. [27]. Equivalently we compute the number 
of ^-induced muon events. 

Neutrino oscillations introduced by NP effects result in an energy dependent 
distortion of the zenith angle distribution of ATM muon events. We quantify this 
effect in IceCube by evaluating the expected angular and E^ m distributions in the 
detector using Eqs. (9) and (10) in conjunction with the fluxes obtained after 
evolution in the Earth for different sets of NP oscillation parameters. 

For illustration we concentrate on oscillations resulting from VLI that lead to 
an oscillation wavelength inversely proportional to the neutrino energy. The results 
can be directly applied to oscillations due to VEP. We show in Fig. 2 the zenith 
angle distributions for muon induced events for different values of the VLI parameter 
5c/ c and maximal mixing £ v u = tt/A for different threshold energy E^ n > -EWcshoid 
normalized to the expectations for pure Am 2 oscillations. The full lines include both 
the ^-induced events (Eq.(9)) and i/ r -induced events (Eq.(10)) while the last ones 
are not included in the dashed curves. We see that for a given value of 5c/ c there is 
a range of energy for which the angular distortion is maximal. Above that energy 
the oscillations average out and result in a constant suppression of the number of 
events. Inclusion of the ^--induced events events leads to an overall increase of the 
event rate but slightly reduces the angular distortion. 

In order to quantify the energy-dependent angular distortion we define the vertical- 
to-horizontal double ratio 
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Fig. 2. Upper panels: Zenith angle distributions for muon induced events for different val- 
ues of the VLI parameter 5c/c and maximal mixing ^ v u = n/A for different threshold energy 
> -^threshold normalized to the expectations for pure Am 2 oscillations . Lower panels: 
The predicted horizontal-to- vertical double ratio in Eq.(10) for different values of <5c/c. The data 
points in the figure show the expected statistical error corresponding to the observation of no NP 
effects in 10 years of IceCube. 



NV U {E fi 



-0.6 < cos 6» < -0.2) 



hor 



N 



no~vli 



-0.6 < cost? < -0.2) 



-1 < cosfl < -0.6) 



(10) 



/' 



— vli 



■1 < cost? < -0.6) 



where by Ej L m ' 1 we denote integration in an energy bin of width 0.2 log 10 ( J B^ m ' 1 ) 
using that IceCube measures energy to 20% in log 10 E for muons. 

In what follows we will use the double ratio in Eq. (10) as the observable to 
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determine the sensitivity of IceCube to NP-induced oscillations. We have chosen 
a double ratio to eliminate uncertainties associated with the overall normalization 
of the ATM fluxes at high energies. It is worth noticing that using this observable 
relies on the fact that the zenith angular dependence of the effective area is well 
understood. 

In Fig. 2 we plot the expected value of this ratio for different values of 8c/ c. As 
mentioned above, IceCube measures energy to 20% in log 10 E for muons. Accord- 
ingly, we have divided the data in 16 E^ m bins: 15 bins between 10 2 and 10 5 GeV 
and one containing all events above 10 GeV. In the figure the full lines include 
both the ^-induced events (Eq.(9)) and ^-induced events (Eq.(10)) while the last 
ones are not included in the dashed curves. As described above, the net result of 
including the z/ r -induced events is a slight decrease of the maximum expected value 
of the double ratio. The data points in the figure show the expected statistical error 
corresponding to the observation of no NP effects in 10 years of IceCube. In order 
to estimate the expected sensitivity we assume that no NP effect is observed and 
define a simple x 2 function including only the statistical errors. 

We show in Fig. 3 the sensitivity limits in the [Sc/c, £ v n] -plane at 90, 95, 99 and 
3 a CL obtained from the condition x 2 (5c/ c, £ v ii) < Xmax (CL, 2dof). We show in 
the figure the results obtained using the RQPM model and the TIG model. The 
difference is about 50% in the strongest bound on 5c/ c. 




Fig. 3. Sensitivity limits in the <5c/c,£ v i; at 90, 95, 99 and 3 <r CL. The hatched area in the upper 
right corner is the present 3<r bound from the analysis of SK data in Ref. [16]. 

The figure illustrates the improvement on the present bounds by more than two 
orders of magnitude even within the context of this very conservative analysis. The 
loss of sensitivity at large 5c/ c is a consequence of the use of a double ratio as an 
observable. Such an observable is insensitive to NP effects if 5c/ c is large enough 
for the oscillations to be always averaged leading only to an overall suppression. 

When data becomes available a more realistic analysis is likely to lead to a further 
improvement of the sensitivity. 
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